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ABSTRACT 

Aims. We investigate the temporal spot evolution of the K-giant component in the RS CVn-type binary system f Andromedae to 
establish its surface differential rotation. 

Methods. Doppler imaging is used to study three slightly overlapping spectroscopic datasets, obtained independently at three different 
observing sites. Each dataset covers one full stellar rotation with good phase coverage, and in total, results in a continuous coverage 
of almost three stellar rotations (P mt = 17.8 d). Therefore, these data are well suited for reconstructing surface temperature maps and 
studying temporal evolution in spot configurations. Surface differential rotation is measured by the means of cross-correlation of all 
the possible image pairs. 

Results. The individual Doppler reconstructions well agree in the revealed spot pattern, recovering numerous low latitude spots 
with temperature contrasts of up to «1000K with respect to the unspotted photosphere, and also an asymmetric polar cap which 
is diminishing with time. Our detailed cross-correlation study consistently indicate solar-type differential rotation with an average 
surface shear a w 0.055, in agreement with former results. 
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1. Introduction 

Stellar activity involves a range of magnetic activity phenomena, 
which are generated and modulated by the dynamo mechanism. 
The differential rotation is of utmost importance in understand- 
ing the magnetic activity, since it is a key element of the stellar 
dynamo, which has a controlling influence over the strength of 
magnetic fields generated, and thus that of the activity itself. The 
differential rotation on stars with convective envelopes puts con- 
straints on the large scale topology of the magnetic field, there- 
fore giving important information on the working of the dynamo 
beneath the surface. Nevertheless, on stars the detection of the 
surface differential rotation still remains a challenging obser- 
vational task, since stellar surfaces, except in rare cases (e.g., 
Gilliland & Dupree [l996l Monnier et al. l2007l . cannot yet be re- 
solved directly. Thus, indirect surface reconstruction technique, 
like Doppler imaging (see, e.g., Vogt et al. I19871 Piskunov et 
al. 1 1990b . is needed. Moreover, since starspots are commonly 
regarded as tracers of photospheric plasma motions, the pattern 
of surface differential rotation can be evolved by investigating 

Send offprint requests to: Zs. Kovari, e-mail: kovari@konkoly.hu 

* Based on observations collected using the Bernard Lyot Telescope 
at Pic du Midi Observatory; the 1.93-m telescope at Haute-Provence 
Observatory, France; and the 8.2-m Kueyen telescope (VLT/UT2) of 
the European Southern Observatory, Chile (Prg. 081.D-0216(A)). 



short-term spot redistribution through reconstructed time-series 
surface images (e.g., Donati & Collier Cameron 1 19971 Weber & 
Strassmeier l20fjn K ovari e t al. 120041 l2007al Barnes et al.2005, 
Hussain et al . 120061 Weber l2007l etc.). 

In this work the differential rotation of the K-giant compo- 
nent in the long -period RS CVn binary f Andromedae (hereafter 
(And) is studied. (And (HD4502) is a single-lined spectro- 
scopic binary with a bright (V = +4.1 from SIMBAD) K-giant 
ellipsoidal variable component (Campbell et al. 1191 11 1, and an 
unseen companion of likely late-G or early-K dwarf (Kovari 
et al. I2007bl hereafter Paper I, see also the references therein). 
Although, its photometric variability is ruled by the ellipticity 
effect, light curve modulation with the rotational cycle of 17.8 
days is also present, suggesting spot activity (cf. Strassmeier 
et al. I1989I I. This feature is in agreement with the strong Can 
H&K emission reported first by Joy & Wilson (1949). The star 
is known also from its Ha, UV, and X-ray activity (cf. Eaton 
1995 Reimers fl980l and Schrijver et al l 19841 respectively). 

The first extended study of the target using photometric mod- 
elling and Doppler imaging was presented in Paper I. Doppler 
maps revealed spots dominating at lower to medium latitudes 
with a typical temperature contrast of ^750 K, as well as a 
weak and variable polar cap. Considering the distorted geom- 
etry, the fundamental system and stellar parameters were speci- 
fied and an investigation for differential rotation was presented. 
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Recently, Korhonen et al. (120101 Paper II from here on) mea- 
sured the angular diameter of the star with the Very Large 
Telescope Interferometer (VLTI/ AMBER) and specified its ab- 
solute dimensions. A new Doppler image was also reconstructed 
from the high quality spectra obtained simultaneously from the 
VLT/UVES instrument. The resulting image, similarly to the 
ones in Paper I, recovered mainly cool spots at lower lati- 
tudes, amongst with the most dominant features concentrating 
at quadrature positions, and also some polar features, but with 
significantly weaker contrast. 

A unique opportunity presented itself during August and 
September 2008, when the star was measured almost con- 
secutively from Pic du Midi Observatory (France), from 
Haute-Provence Observatory (France), and from the European 
Southern Observatory (Chile, see Paper II), providing three in- 
dependent (slightly overlapping) spectroscopic datasets, each of 
them suitable to reconstruct one single Doppler image using 
four favorable mapping lines. In Sect. [2] a brief review of the 
three datasets is given, then, in Sect. [3] the results from Doppler 
imaging are presented. Afterwards, in Sect. |4]the Doppler maps 
are compared with each other by using our sophisticated cross- 
correlation technique to search for the surface differential rota- 
tion pattern. Finally, the results are summarized and discussed in 
Sect.|5] 

2. Spectroscopic observations 

The spectroscopic data were collected at three different observ- 
ing sites independently, between 13 August and 1 October, 2008. 
Each dataset covers roughly one stellar rotation (16-18 days) 
and the consecutive sets overlap one another by 1-2 days. In 
all the three cases, the covered wavelengths include the 6400 A 
region commonly used for Doppler imaging. This contains, 
amongst others, the mapping lines Fe 1-6411, Fe 1-6421, Fei- 
6430 and Ca 1-6439. The summary of the observations is given 
in TableQ] 

All the spectra are phased using the same ephemeris as in 
Paper I, fitting the zero phase to the conjunction with the sec- 
ondary in front: 

HJD = 2449 997.223 ±0.017+ 17.769426 ± 0.000040 x E. (1) 

2.1. NARVAL spectra 

The first dataset in the time-series with altogether 48 spectra was 
obtained during 10 nights between 13-31 August at Pic du Midi 
Observatory (France) using the 2-m Bernard Lyot Telescope, 
equipped with the NARVAL spectrograph (Auriere 2003). The 
spectrograph covers the wavelength range between 3700 A and 
10000 A with a resolving power (A/AA) of 60000. Each observa- 
tion consists of four 30 s long exposures obtained immediately 
after each other. These exposures were combined into one ob- 
servation to increase the signal-to-noise ratio (S/N). This way 
10 individual observations (one per night) from different rota- 
tional phases were obtained. The information in Table Q] is for 
the combined observations. The reference HJD of the Doppler 
map calculated from the NARVAL spectra in Sect. [3] (i.e., the 
mean HJD) is 2454702.088. 

2.2. SOPHIE spectra 

The second set consisting of 10 spectra were collected between 
30 August and 15 September, 2008 with the SOPHIE echelle 



Table 1. Observational records of our spectroscopic data. 



HJD 


phase 


S/N 


date 


instrument 


n b 


t c 


2454000+ 






(2008) 






Tsl 


691.5464 


0.178 


446 


13 Aug 


NARVAL 


4 


30 


696.5723 


0.461 


258 


18 Aug 


NARVAL 


4 


30 


699.6046 


0.632 


242 


21 Aug 


NARVAL 


4 


30 


701.5989 


0.744 


679 


23 Aug 


NARVAL 


4 


30 


702.5645 


0.798 


384 


24 Aug 


NARVAL 


4 


30 


703.5135 


0.852 


270 


25 Aug 


NARVAL 


4 


30 


704.5673 


0.911 


296 


26 Aug 


NARVAL 


4 


30 


705.5435 


0.966 


329 


27 Aug 


NARVAL 


4 


30 


708.6240 


0.139 


421 


30 Aug 


SOPHIE 


1 


120 


708.6278 


0.140 


371 


30 Aug 


NARVAL 


4 


30 


709.5420 


0.191 


445 


31 Aug 


SOPHIE 


1 


120 


709.6025 


0.194 


291 


31 Aug 


NARVAL 


4 


30 


710.5015 


0.245 


335 


01 Sep 


SOPHIE 


1 


120 


711.6389 


0.309 


441 


02 Sep 


SOPHIE 


1 


120 


714.6815 


0.480 


261 


05 Sep 


SOPHIE 


1 


120 


715.6242 


0.533 


198 


06 Sep 


SOPHIE 


1 


950 


717.4957 


0.639 


269 


08 Sep 


SOPHIE 


1 


120 


718.5432 


0.698 


441 


09 Sep 


SOPHIE 


1 


120 


720.5707 


0.812 


237 


11 Sep 


SOPHIE 


1 


166 


722.8167 


0.940 


866 


13 Sep 


UVES 


3 


8 


724.5445 


0.035 


258 


15 Sep 


SOPHIE 


1 


120 


11 A l(\Af\ 


U.U'+O 


ash 




U V Co 




£ 

o 


726.7101 


0.159 


586 


17 Sep 


UVES 


3 


8 


727.6870 


0.214 


628 


18 Sep 


UVES 


3 


8 


728.7387 


0.273 


762 


19 Sep 


UVES 


3 


8 


730.7203 


0.384 


693 


21 Sep 


UVES 


3 


8 


732.7287 


0.497 


914 


23 Sep 


UVES 


3 


8 


734.6517 


0.606 


785 


25 Sep 


UVES 


3 


8 


736.6936 


0.721 


612 


27 Sep 


UVES 


3 


8 


740.7533 


0.949 


802 


01 Oct 


UVES 


3 


8 



Notes. (<!) Phases computed using Eq.Q] (,,) Number of exposures com- 
bined into a final spectrum. <c) Lengths of the single exposures. 



spectrograph (Perruchot et al. 120081 1 at the 1.93-m telescope of 
Haute-Provence Observatory (France). The spectrograph has a 
wavelength range of 3872-6943 A. The observations were car- 
ried out in the high efficiency mode (HE) using 100// fibre, 
which gives A/AA of 40000. The exposure time was usually 
120 s, but in two cases longer times were used: on September 
6 the exposure time was 950 s and on September 11 166 s. The 
mean HJD of the SOPHIE observations is 2454715.177. 



2.3. UVES spectra 

The third set of observations were carried out at the European 
Southern Observatory (Ch ile) wi th UVES (UV- Visual Echelle 
Spectrograph, Dekker et al. 2000) spectrograph mounted on 8.2- 
m Kueyen VLT unit telescope. The data were collected during 
ten nights between 13 September and 1 October, 2008. In the 
observations the red arm standard wavelength setting of 600 nm 
was used with image sheer #3. This configuration gives a re- 
solving power (A/AA) of 1 10000 and a wavelength coverage of 
5000-7000 A. For every observation the three 8 s exposures ob- 
tained immediately after each other were combined into a sin- 
gle, very high S/N spectrum. For further details see Paper II. 
Computing the average HJD we get 2454730.621. 
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3. Doppler imaging with TempMap £ 

All of the three datasets cover roughly one stellar rotation, and 
thus allow to reconstruct one Doppler image for each set. For 
the surface image reconstruction we started from the Doppler 
imaging code TempMap (Rice et al. [T989) which, afterwards, 
was extended by us (TempMap,;, Paper I) to take into account 
the non-spherical shape of stars in close binary systems. In our 
approximation, the real Roche-geometry is simplified to a rota- 
tion ellipsoid, that is elongated towards the other component. We 
note here that TempMap £ is the only Doppler imaging code avail- 
able so far, that can use non-spherical surface geometry. We also 
note, that neglecting the ellipticity in the Doppler reconstruction 
would yield a significant rise in the goodness-of-fit value of the 
line profile inversion (see Papers I-II). Eventually, the astrophys- 
ical input data of f And, including the most feasible distortion 
parameter, are adopted from Paper I (Table 2 therein). 

Doppler imaging was performed for a sample of well known 
mapping lines, that usually produces the most reliable outcome. 
In practice, we use only those lines that are jointly covered by 
all the three datasets: Fe 1-6411, Fe 1-6421, Fe 1-6430 and Cai- 
6439. The resulting Doppler images are plotted in Fig.Q] Albeit 
a higher resolving power would allow for a larger amount of 
grid points to be used in the map, for the comparison reasons 
we apply the same grid size of 5° x 5° in all the maps. The in- 
dividual reconstructions well agree in the revealed spot pattern, 
i.e., numerous low latitude spots are recovered with tempera- 
ture contrasts of up to »1000K with respect to the undisturbed 
photosphere of 4600 K. In addition, an asymmetric cool polar 
feature is also present at around <p « 0.75, however, switching 
from the NARVAL map to the SOPHIE and UVES maps, i.e., 
from August to September 2008 the contrast became weaker and 
the spot concentration shifted backwards in phase. We refer also 
to Paper II, wherein, using the UVES observations, a thorough 
Doppler imaging study was carried out. 



4. Cross-correlation study 

Surface differential rotation (DR) of a spotted star can be 
measured by cross-correlation of subsequent Doppler images 
(Donati & Collier Cameron |1997| l. Nevertheless, even when hav- 
ing reliable Doppler images, short-term spot changes (e.g., vivid 
interaction between neighbouring spots, emerging a new spot 
or disintegrating an old one, etc.) can easily mask the trace 
of the surface DR in the cross-correlation function (hereafter 
ccf) map. To emphasize the DR pattern, and to diminish differ- 
ences originating from individual spot motion and evolution, we 
use our newly developed cross-correlation technique ACCORD 
('Average Cross-CORrelation of time-series Doppler images', 
see e.g., in Kovari et al. 2004, 2007a, Paper I). The basic idea 
of the method is, that averaging the more ccf maps, the better re- 
sult can be achieved by boosting such cross-correlation features 
that are jointly present, like the DR pattern itself. 

From the Doppler maps shown in Fig. Q] three cross- 
correlation pairs can be formed: NARVAL/SOPHIE (N/S), 
NARVAL/UVES (N/U), and SOPHIE/UVES (S/U). Those ccf 
maps then can be averaged. However, a simple averaging will not 
definitely boost the DR pattern, since the shape of the correlation 
pattern depends on the time gap between the paired Doppler im- 
ages, i.e., the larger the time gap, the larger the elongation of 
the individual 1-D longitudinal cross-correlation functions par- 
allel to the equator (in practice, these are longitude strips in bins 
of 5°, see, e.g., Kovari et al. I2004I I. Thus, to avoid fuzziness, 
a linear renormalization must be done before averaging. When 
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Fig. 1. From top to bottom: the combined (Fe 1-641 1+Fei- 
6421+ Fe i-6430+Ca 1-6439) Doppler maps in time series for the 
NARVAL, SOPHIE, and UVES data, respectively. Arrows un- 
derneath indicate the observational phases. 



the final average ccf map is assembled, the best-correlating pat- 
tern can be fitted by an assumed form of the DR law. The result 
of applying ACCORD to the three Doppler images is plotted in 
Fig. [2] For the latitude (J3) dependent rotation law we assume the 
usual quadratic form of Q.(J3) - Cl eq (1 - asm 2 fi), where Q. eq is 
the equatorial angular velocity and a = (£2 eq - £2 po i e )/Q eq is the 
dimensionless surface shear parameter. The best fit yields solar- 
type equatorial acceleration: 

Q08) = 20.763 (1 - 0.053 sin 2 /?) [°/d]. (2) 



5. Discussions and conclusion 

The outcome of our new cross-correlation study seems to be 
fairly consistent with the result in Paper I, where a of 0.049 + 
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Fig. 2. Applying the ACCORD technique for the three available 
Doppler images plotted in Fig. Q] Better correlation on the re- 
sulting ccf map is represented by darker shade. The average 
longitudinal cross-correlation functions in bins of 5° are fitted 
by Gaussian curves (Gaussian peaks are indicated by dots, the 
FWHMs by horizontal lines). The continuous line represents the 
best fit quadratic DR law. 

0.003 was derived. To test the reliability of our method and re- 
sult we repeat the correlation study, but instead of using the com- 
bined Doppler images in Fig.[T| we start from individual recon- 
structions, i.e., we use Doppler images that are recovered indi- 
vidually for a single mapping line. This way we have four maps 
(Fe i-64 1 1 , Fe 1-642 1 , Fe 1-6430 and Ca 1-6439) for each datasets 
(N, S, U), separately, thus 12 cross-correlations can be achieved 
and finally averaged. We expect that this way the signal-to-noise 
of the DR pattern in the ultimate average ccf map will improve. 
The result with the fitted DR pattern is shown in Fig. [3] yielding 
nearly the same DR law as of Eq. [2] 

£206) = 20.689 (1 - 0.055 sin 2 /?) [°/d], (3) 

nevertheless, the relative height of the fitted 'ridge' in the ccf 
landscape (i.e., the maxima of the gaussian peaks) has signifi- 
cantly (^40% on average) risen, compared to the one in Fig. [2] 
We can also estimate the errors of the derived parameters given 
in Eq. [3] by carrying out the whole process, but with using 
only three mapping lines from the available four, i.e., we per- 
form four cross-correlation studies separately and estimate the 
error bars from the 1-cr of the obtained DR parameters, getting 
e(Q eq ) = +0.056 °/d and s(a) = +0.0022. We consider this result 
as a proof of the reliability and robustness of the ACCORD tech- 
nique. A comparative summary of the derived DR parameters is 
given in Table [2] 

Theoretical calculations predict that across the H-R diagram 
DR becomes stronger with increasing rotation period and evolu- 
tionary status (Kitchatinov & Riidiger 1999). Still, £ And shows 

Table 2. Comparing the DR parameters 



[°/d] 


a 


LT 7 
[d] 


note 


19.027 


0.049 


381 


from Paper I 


20.763 


0.053 


332 


using combined images 


20.689 


0.055 


315 


using separate reconstructions 



Notes. <a) Time needed by the equator to lap the pole by one full period. 
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Fig. 3. Applying the ACCORD technique for the individual 
Doppler reconstructions using Fe 1-6411, Fe 1-6421, Fe 1-6430 
and Ca 1-6439 mapping lines, separately. Otherwise as in Fig. [2] 

relatively weak surface differential rotation. This could be ex- 
plained by large starspots not following the stellar surface flows, 
as they are anchored at different depth in the convection zone. 
Recent investigation by Korhonen & Elstner (1201 II ) using dy- 
namo models shows that the much weaker DR than what is 
used in the model is recovered when having only the large-scale 
dynamo field. The real input DR law can only be recovered 
when the model also includes small-scale fields. The question 
is whether or not large starspots can really be caused by small- 
scale magnetic fields. 

According to Kitchatinov & Riidiger (2004) giant stars in 
close binary system, e.g., RS CVn-type binaries, can also show 
anti-solar differential rotation, where the equator rotates slower 
than the poles. Such a DR law has been detected, for example, 
in o-Gem (Kovari et al. l2007al . Kitchatinov & Riidiger (12004b 
attribute the anti-solar differential rotation to strong meridional 
flow, which could be caused by barocline driving due to large 
thermal spots or by field forcing of a close companion. In accor- 
dance with this, a poleward meridional flow has been reported 
for cr Gem with an average velocity of 350 km s (Kovari et al. 
2007a). On the other hand, no sign of such a meridional flow is 
seen on £ And (Paper I). This observed difference between £ And 
and cr Gem can be rooted in the different system geometry and 
dynamics (Kovari et al. 2012). Distortion and tidal forces could 
account for the converse DR in these two cool giants. 
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